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а 要 : 光 干 涉 要 求实 现 各 光路 共 相 位 登 加 , 对 于 宽 波 带 的 白光 干涉 要 求 则 更 加 苛刻 。 因此 ， 
要 实现 稳定 的 白光 光 王 涉 需 要 对 光 程 进行 实时 探测 和 控制 。 本 文 对 基于 传 里 叶 变换 的 色散 条 纹 
法 进行 了 理论 推导 、 算 法 设计 ， 并 利用 实验 室 建立 的 Fizeau 型 干涉 实验 装置 开展 开 、 闭 环 实验 
研究 。 在 开 环 实验 中 ， 采 集 不 同 光 程 差 的 多 幅 色 散 条 纹 图 像 ， 研 究 了 图 像 频谱 次 峰 偏 移 量 和 光 
程 差 的 关系 ， 结 果 显 示 ， CEE а ида ер 在 闭环 实验 中 ， 利 用 光 
实时 计算 结果 和 光 程 补偿 机 构 ， 对 光 程 差 进行 闭环 控制 ， 结 果 显 示 ， 在 外 加 扰动 的 情况 

系统 能 够 始终 保持 在 初始 的 干涉 状态 。 

Вал 光 程 探测 ; 傅 里 叶 变 换 ; 色散 条 纹 法 ; 控制 
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光学 综合 孔径 技术 是 将 多 个 小 口径 的 光学 元 件 或 系统 以 某 种 构 型 进行 精确 排列 ， 使 通过 各 
子 孔 径 的 光束 在 焦 平 面 上 实现 光 场 的 相干 阁 加 ， 经 图 像 复 原 技术 处 理 后 ， 达 到 与 单一 大 口径 系 
统 相 当 的 衍射 极限 分 辨 率 的 技术 。 它 可 以 突破 单 望远镜 口径 衍射 极限 的 限制 。 

光学 综合 孔径 望远镜 按照 拼接 方式 的 不 同 , 分 为 拼接 主 镜 和 拼接 子 望远镜 两 种 形式 。 目 前 ， 
世界 各 国 已 建立 了 多 台 拼 接 子 望远镜 形式 的 光学 综合 孔径 望远镜 。2002 F, 美国 床 省 理工 学 
院 建立 了 白光 自 适 应 GOLAY-3 型 ARGOS 望远镜 "， 它 由 3 个 孔径 为 21cm 的 子 望 远 镜 组 成 ， 角 分 
HIK 0. 35"， 光 谱 范 围 400-700nm， 视 场 为 3'X3'， 等 效 口 径 0. 62m， 信 噪 比 100。 美 国 洛克 希 
德 马丁 公司 联合 NASA, JL 实验 室 、 加 州 理工 学 院 、 加 利 福 尼 亚 大 学 和 罗切斯特 大 学 ， 建 
立 了 空间 多 仪器 口径 分 布 式 传 感 系统 MIDAS) ”, 该 系统 能 够 实现 大 视 场 、 衍 射 极限 成 像 ， 
可 以 在 5000km 高 的 轨道 上 进行 对 地 观察 ， 地 面 分 辨 率 小 于 1m, 也 可 以 在 100km 高 的 轨道 上 进 

行 高 分 辨 率 成 像 ， 地 面 分 辨 率 达 到 2cm。 

光学 综合 孔径 望远镜 在 进行 干涉 成 像 时 ， 必 须 实 现 光 束 的 相干 登 加 ， 才 能 提高 图 像 的 空间 
分 辩 率 。 对 于 和 白光 宽带 成 像 系 统 ， 由 于 光源 的 相干 长 度 非 常 小 ， 相 当 于 要 求 不 同 路 径 的 光束 实 
现 零 等 光 程 干涉 。 因 此 ， 对 光 程 差 进 行 精确 探测 ， 就 成 为 了 实现 光 王 涉 的 先决 条 件 。 
目前 ， 国 内 外 在 光 程 探测 领域 已 开展 了 大 量 的 研究 。 出 瞳 面 检测 : 四 棱锥 波 前 检测 法 通过 
探测 出 瞳 处 4 个 像 点 的 能 量 差 异 计算 光 程 差 ”。 离 焦 面 检测 : 相位 恢复 法 通过 焦 面 和 离 焦 面 的 
一 对 光 强 分 布 数据 反 解 出 各 自 的 相位 分 布 ”。 焦 面 检测 法 : (D) 罕 带 一 哈 特 曼 夏 克 法 和 宽带 一 哈 
特 曼 夏 克 法 通过 检测 采样 孔径 衍射 图 样 和 “样板 图 样 ” 的 相关 系数 ， 计 算 拼 接 子 镜 间 的 光 程 差 
本。(2) 色 散 条 纹 法 通过 检测 所 提取 的 相 邻 孔径 之 间 的 干涉 条 纹 中 横向 光 强 信号 ， 对 其 进行 非 线 
性 最 小 二 乘 拟 合 ， 求 解 光 程 差 "。 

ЛАН 8 达到 光 程 差 探测 范围 大 于 +50pm, 探测 精度 小 于 20nm 的 技术 指标 , 并 利用 国家 天 广 
台 空 间 天 文 技术 实验 室 所 建立 的 光学 综合 孔径 原理 样机 系统 ， 本 文 对 基于 傅 里 时 变换 的 色散 条 
纹 法 进行 了 理论 推导 和 算法 设计 ， 并 开展 了 开 环 和 闭环 实验 。 在 开 环 实验 中 ， 采 集 不 同 光 程 差 
的 色散 条 纹 图 像 ， 通 过 图 像 频谱 次 峰 偏 移 量 解 算 光 程 差 ， 并 与 设 定 的 光 程 差 进行 比 对 ， 验 证 方 
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程 差 变化 的 校正 能 


1 理论 基础 


法 的 线性 关系 ;在 闭环 实验 9 


FPF， 对 人 为 引入 的 光 程 差 进 行 实 时 计算 和 闭环 控制 ， 观 察 系统 对 光 


基于 色散 条 纹 法 的 焦 平 面 光 强 分 布 表达 式 为 ""; 
ee | (1) 
ns 
为 色散 方向 ，C 为 线 色散 率 ，g,(y) 为 初 相位 。 由 于 波 数 的 表达 式 K(x) = 代 将 
си PETF 


该 表达 式 进 行 1 阶 泰勒 近似 展开 ， 得 至 


Mk=. FK) В ЗЕ 


0 


近似 展开 式 带 入 (1)， 可 得 : 
ky C, 
I(x, у)< 1 + 7-сов (ko — а х)-ф+ P(Y) (2) 
对 式 (2) 进 行 二 维 傅 里 叶 变 换 ， 得 到 : 
—j2m + уй, 
Lo f) = TG, ye F (E7 УЛ dydy 3) 


由 于 表达 式 (2) 中 的 余弦 项 可 以 月 


св, БО 


Ilfo f) я 2 90) + пре “9 + 


СА. 


+rhy e” :д(.- 


рани лае 


печки 
表达 式 为 : 


k 


日 欧 拉 公式 表达 ， 因 


Ох). sao) 


因此 7(x, y) 中 余弦 项 可 表示 为 : 


ве 


Ох). Hao) дема y). rmt 
е 


+— 4 
2 2 (4) 


руни ar rA) 


Ва). E 


k и: -А) o 


] 式 (4) 来 表示 (2) 中 的 余弦 项 ， 得 到 了 (x, y) 的 二 维 频谱 表达 式 : 


C ша 


9, f,—A) 


(5) 


E Tx, y) 的 频谱 中 有 三 个 峰 , 其 中 主峰 位 于 (0,0)， 


因此 , 其 中 一 个 次 峰 相 对 于 y 轴 的 偏 移 量 dx 
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dx = 一 人 (6) 


其 中 Cu 和 均 为 常数 , 因此 频谱 次 峰 相 对 于 y 轴 的 偏 移 量 必 与 光 程 差 儿 之 间 呈 现 出 良好 的 线性 
关系 ， 故 可 以 用 基于 傅 里 叶 变 换 的 方法 对 色散 条 纹 图 像 进行 处 理 ， 从 而 求 得 光 程 差 值 。 


2 光 程 探测 算法 流程 设计 


基于 色散 条 纹 图 像 的 光 程 探 测算 法 流程 如 图 1 所 示 ”。 首 先 ， 读 取 某 一 光 程 差 下 的 干涉 条 纹 
图 像 ， 然 后 对 其 进行 二 维 全 里 叶 变换 ， 并 求 取 幅 值 。 得 到 频谱 图 像 后 ， 首 先 将 主峰 所 在 的 区 域 
全 部 置 零 ， 然 后 求 取 最 大 幅 值 所 在 的 位 置 ， 此 时 所 求 得 的 就 是 两 个 次 峰 所 在 的 位 置 。 选 取 其 中 
一 个 次 峰 ， 并 确定 该 次 峰 相对 于 y 轴 的 偏 移 量 dx ， 以 该 次 峰 所 在 位 置 为 中 心 ， 进 行 亚 像素 精度 
Ма, ВЕР (Бу), ИНА ЕО у), MERA”: 


05 ТЕСх-1,у)-Е(х+1,у)| 
[F(x-1, y)+ F(x+1, y)—2F (x, y)] H 
在 此 基础 上 ， 乘 以 一 个 比例 系数 就 可 以 求 得 光 程 差 值 ， 比 例 系 数 可 以 通过 计算 得 到 也 可 由 实验 
进行 标定 ， 至 此 ， 对 于 一 幅 干涉 条 纹 图 像 的 光 程 差 求解 过 程 结束 。 
C 开始 ) 


Iar o 
FFP AA ва ва 


戏 [二 涉 条 纹 图 像 进行 二 维 博 里 时 变换 | 


ах = 4х+ 


| REA | 


| 主峰 所 在 范围 置 零 | 


| 以 其 中 一 个 次 峰 所 在 的 行 和 列 | 


| 塘 访 次 峰 相 对 于 y 轴 的 整数 倍 次 峰 位 移 | 


| чи (Фр зе жу ни ол Фр | 


C ÆR ) 
图 1 基于 一 幅 图 像 的 光 程 探测 算法 流程 图 


Fig.1 Piston error detection algorithm flow based on one image 
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国家 天 文 台 建 立 的 Fizeau 型 光学 综合 孔径 实验 装置 由 光源 系统 、 子 望远镜 系统 、 指 向 探测 
和 补偿 系统 、 光 程 探 测 和 补偿 系统 、 合 光 成 像 及 处 理 系统 5 部 分 组 成 ， 见 图 2”。 其 电子 学 系 
统 由 图 像 分 发 板 、3 块 计算 控制 板 、 数 模 输 出 板 、 模 拟 电源 板 、 数 字 电 源 板 和 接口 底板 共 8 块 


板 卡 组 成 ， 见 图 3。 在 开 环 实验 中 ， 采 集 1 


色散 条 纹 的 系列 图 像 ， 利 用 MATLAB 求解 光 程 差 。 在 


闭环 实验 中 ， 利 用 电子 学 系统 ， 基 于 Visual Dsp++5.0 开发 TS201 的 程序 ， 进 行 光 程 差 的 实时 探 


测 及 控制 。 


在 实验 中 ， 光 源 系 统 均 采 用 NKT Photonics 公司 的 SuperK EXTREME 系列 高 功率 超 连 续 白 


光 光 源 ， 其 滤 光 器 输出 光谱 范围 


为 400nm 至 700nm， 


有 具有 波段 中 心 、 波 段 宽度 、 输 出 功率 可 调 


的 特点 。 光 源 望 远 镜 是 典型 的 卡 塞 格林 望远镜 ， 具 有 500mm 通 光 孔径 ，20m 焦距 ， 大 于 2arcmin 


的 视 场 角 ，F#/40。 子 望远镜 系统 上 


h3 路 子 望 远 镜 组 成 ， 采 用 反射 式 无 焦 系 统 ， 放 大 率 为 5， 基 


线 长 度 200mm， 入 射 光 束 口 径 100mm， 出 射 光束 口径 20mm， 系 统 组 合 焦距 4500mm。 为 保证 经 子 
望远镜 lH, 2#, 3# 的 三 路 光束 可 进行 相干 成 像 ， 光 程 探 测 模块 内 设置 2 个 光 程 探测 望远镜 ， 以 子 
望远镜 1# 为 基准 ， 分 别 对 1#，2# 和 和 1#，3# 进 行 光 程 探测 并 解 算 出 光 程 差 值 ， 最 终 达 到 控制 3 
路 光束 间 光 程 差 的 目的 。 其 中 1#，2# 和 1#，3# 间 的 两 路 光束 经 光 程 探测 模块 内 的 Amici 色散 楼 
机 以 及 光 程 探测 望远镜 后 ， 在 焦 面 上 形成 色散 干涉 条 弘 。 光 程 补偿 模块 分 为 粗 调 平 台 与 精 调 平 
， 利 用 光 程 补偿 组 件 中 的 粗 调 平 台 〔PI 公司 的 M-122 微 位 移 平 


台 ， 如 图 4 所 示 。 在 开 环 实验 中 


п) 调整 运动 反射 镜 的 位 置 改变 光学 延迟 线 的 大 小 ， 从 而 得 到 不 同 光 程 差 下 的 干涉 条 纹 图 像 。 


M-122 微 位 移 平 台 以 线性 编码 器 作为 集成 传感器 , 行程 5mm, 分 辨 率 0.1pm, 最 小 位 移 0.2hm， 


最 大 速度 20mm/s。 在 闭环 实验 


110mm， 焦 距 1680mm。 在 完成 倾斜 误差 与 光 程 差 的 校正 后 ， 该 系统 将 3 路 光束 进行 相干 成 像 。 


图 2 原理 样机 系统 


FP， 利用 光 程 补偿 组 件 中 的 精 调 平台 (PI 公司 的 P-753 促 动 器 )， 
对 光 程 差 进行 精 补偿 。P-753 促 动 嚣 带 有 电容 式 集成 传感器 ， 行 程 12pm， 分辨 率 0.05nm， 重 复 
精度 土 Inm。 合 光 成 像 系统 是 由 1 个 合 光 成 像 望 远 镜 组 成 ， 该 望远镜 为 卡 赛 格 林 形 式 ， 口 径 


Fig.2 Principle prototype System 


4 开 环 实验 


图 4 光 程 补偿 组 件 粗 调 平 


Fig.4 Piston compensation elements coarse and fine adjustment platform 


图 3 电子 学 系统 


Fig.3 Electronic system 


台 和 精 调 平台 
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出 的 光谱 范围 为 550-650nm、 中 心 波长 为 600nm。 通 


过 移动 M-122 平台 ， 在 -90hm 至 +90pm 的 光 程 差 范围 内 ， 以 1pm 为 间隔 ， 采 集 两 束 光 干涉 所 


形成 的 181 幅 色 散 干 涉 条 纹 图 像 。 图 5、 


6 中， 从 上 至 下 分 别 为 -60hm、0hm、+60hm 光 程 差 


下 的 色散 干涉 条 纹 图 像 及 其 频谱 。 当 光 程 差 为 零 时 ， 干 涉 条 纹 未 发 生 倾 斜 。 当 光 程 差 为 正 、 负 


时 ， 干 涉 条 纹 分 别 斜 向 下 、 斜 向 上 倾斜 。 上 


日 于 干涉 条 纹 的 倾斜 线 与 两 个 次 峰之 间 的 连接 线 互相 


垂直 ”, 因此 当 光 程 差 为 零 时 ， 两 个 次 峰 对 称 分 布 于 y 轴 正 、 负 半 轴 。 当 光 程 差 为 正 、 负 时 ， 两 
个 次 峰 对 称 分 布 于 一 、 三 象限 和 二 、 四 象限 。 


图 5 不 同 光 程 差 的 干涉 条 纹 图 像 


Fig.5 Interference fringe images with different piston errors 


图 6 与 干涉 条 纹 图 像 所 对 应 的 频谱 图 


Fig.6 Spectrograms related to interference fringe images 


对 181 幅 图 像 求 得 的 偏 移 量 4 及 其 拟 合 结果 如 图 7 所 示 。 其 中 横 坐标 为 光 程 差 设 定 值 ， 
单位 hm， 纵 坐标 为 所 选取 的 位 于 x 轴 上 方 的 次 峰 相 对 于 y 轴 的 偏 移 量 dx' ， 单 位 为 pixel。 红色 


圆 点 为 与 181 个 光 程 差 所 对 应 的 偏 移 量 ах 


'， 蓝 色 直 线 为 181 个 偏 移 量 的 线性 拟 合 结果 ， 拟 合 


式 为 dx'==0.2536p+0.9104 ， 因 此 可 以 通过 该 表达 式 计算 出 181 个 光 程 差 值 。 图 8 为 实际 求 
解 光 程 差 与 设 定 光 程 差 的 对 比 图 ， 其 中 横 坐 标 为 设 定 光 程 差 , 单位 hm， 纵 坐标 为 实际 解 算 出 的 


光 程 差 值 ， 单 位 为 um。 绿色 圆 点 为 实际 求解 出 181 个 光 程 差 的 分 布 ， 蓝 色 直 线 为 设 定 的 光 程 差 
值 。 所 求解 的 部 分 偏 移 量 А" 及 实际 光 程 差 $ 如 表 1、2 所 示 。 


om 一 -mo 一 


一 Peak displacement with respect to the Piston error 
2 +The linear fitting of the tested 
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图 7 偏 移 量 求解 及 线性 拟 合 结果 


Fig.7 Peak displacements and linear fitting results 


100; 


| + Real piston error 
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图 8 设 定 光 程 差 与 实际 光 程 差 对 比 图 


Fig.8 Set and real piston errors comparison 
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表 1 次 峰 相 对 于 y 轴 偏 移 量 


Table 1 Peak displacements related to y-axis 


-40 -20 0 20 40 
-9.1 -4.0 1.0 6.0 11.0 
表 2 光 程 差 解 算 结果 
Table 2 Piston error calculation results 


-40 -20 0 20 40 


定 光 程 差 值 ， 儿 为 实际 解 算 光 程 差 值 ， 测 量 误差 & 的 表达 式 为 ; 


5 闭环 实验 


с-ф-ф 


在 开 环 实验 中 ， 计 算得 到 181 个 光 程 差 的 均 方 根 误差 E， 为 1.564hm。 


rms 


60 


59.0 


到 -90hm 至 +90hm 范围 内 的 光 程 差 计 算 结 后 ， 我 们 计算 了 均 方 根 误差 值 。 设 纺 为 设 


(8) 


光学 综合 孔径 系统 光 程 探测 及 闭环 控制 框图 如 图 9 押 示 , 其 中 0 为 设 定 的 光 程 差 ， 0 为 
实测 的 光 程 差 。 利 用 基于 传 里 叶 变换 的 光 程 计算 方法 ， 通 过 数字 信号 处 理 器 完成 光 程 差 的 实时 
计算 ， 驱 动 执行 机 构 ， 即 P-753， 通 过 改变 镜面 的 位 置 ， 从 而 对 光 程 差 进 行 补偿 控制 。 控 制 方法 
采用 离散 比例 -积分 -微分 控制 ID)， 其 中 比例 控制 的 动作 与 侦 差 大 小 成 正比 ， 其 作用 是 以 茶 一 


> 


比例 对 稳 态 误差 进 
羊 时 刻 产 和 
有 超前 调节 的 作 上 月 


行 消除 ， 积 分 控制 的 动作 与 偏差 对 时 间 的 积分 成 正比 ， 其 作 ) 


是 消除 每 次 采 
的 余 差 ， 微 分 控制 的 动作 与 偏差 的 变化 速度 成 正比 ， 其 效果 是 阻止 被 调 参数 的 变化 ， 
日 。 增 量 式 PID 表 达 式 为 ; 


u(k)=u(k-1)+ K [e(k)-e(k-1)]+ Ке() + K;le(k)-2e(k -1)+e(k -2)]= 
ulk —-1)+qelk)+qelk-1)+q,elk-2) 


(9) 


其 中 go =K,+K,+K，，4=-K, 2K，，qs =K， Жеищю HA k URRA, ка K K 


分 别 代表 PID 控制 器 的 比例 系数 ， 积 分 系数 和 微分 系数 。 


图 9 光学 综合 孔径 光 程 探测 及 闭环 控制 框图 


Fig 9. Optical synthetic aperture piston error detection and close-loop control diagram 


在 闭环 实验 中 ， 我 们 同样 使 用 NKT 光 源 ， 中 心 波 长 设 为 600nm， 光 谱 范 围 为 550-650nm。 设 
置 K,=0 , K=0.04 ,K=0。 图 10、 图 11、 图 12 分 别 为 调节 到 零 光 程 、2.5hm 光 程 差 以 及 5hm 


光 程 差 时 的 开 、 闭 环 实验 结果 。 当 系统 处 于 开 环 状态 时 ， 干涉 条 纹 图 像 随 着 光 程 差 的 增加 而 发 
生 了 偏 移 。 当 系统 处 于 闭环 状态 时 ,不 管 光 程 差 如 何 变 化 , 干涉 条 纹 始 终 稳定 在 零 光 程 的 位 置 ， 


= 


表明 光 程 差 的 闭环 控制 起 到 了 民 好 的 效 来。 
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图 10 零 光 程 开 环 实验 〈 左 ) 与 闭环 实验 〈 右 ) 对 比 图 
Fig 10. Oum piston error open-loop (ей) and close-loop (right) experiments comparison 
Би” 
Оба 
图 11 2.5pm 光 程 差 开 环 实验 〈 左 ) 与 闭环 实验 ( 右 ) 对 比 图 
Fig 11. 2.5um piston error open-loop (left) and close-loop (right) experiments comparison 
-o Da 
US 


15.5820 
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图 12 Shm 光 程 差 开 环 实验 〈 左 ) 与 闭环 实验 ( 右 ) 对 比 图 


Fig 12. Shm piston error open-loop (left) and close-loop (right) experiments comparison 


6 总 结 


7 ŽA 


本 文 主要 介绍 了 我 们 开展 基于 传 里 叶 变换 的 色散 条 纹 法 进行 光 程 探测 及 开 、 闭 环 实验 结果 
的 情况 。 首 先 ， 我 们 进行 了 理论 推导 ， 得 到 次 峰 相 对 于 y 轴 侦 移 量 与 光 程 差 之 间 的 线性 关系 。 
然后 ， 进 行 了 算法 流程 设计 。 最 后 ， 利 用 实验 室 的 光学 综合 孔径 实验 装置 ， 进 行 了 开 、 闭 环 实 
验 。 在 开 环 实验 中 ， 采 集 多 幅 不 同 光 程 差 的 色散 条 纹 图 像 ， 对 图 像 频 谱 次 峰 偏 移 量 和 光 程 差 的 
关系 展开 研究 ， 结 果 显 示 ， 偏 移 量 与 光 程 差 之 间 呈 现 出 良好 的 线性 关系 ;在 闭环 实验 中 ,利用 
光 程 差 的 实时 计算 结果 和 光 程 补偿 机 构 ， 对 光 程 差 进行 闭环 控制 ， 结 果 显 示 ， 在 外 加 扰动 的 情 
况 下 ， 系 统 能 够 始终 保持 在 初始 的 干涉 状态 。 
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Abstract: For optical interference, to realize the co-phasing effect is important, especially 
rigorous for broadband white light interference. In order to realize stable white light interference, we 
need to make real-time detection and compensation control for the piston error. Through theoretical 
derivation, algorithm design for dispersed fringe sensing based on Fourier transform and open-loop, 
close-loop research with the Fizeau interference experiment facilities in our lab. In the open-loop 
experiment, we collect many dispersed fringe images with different piston errors and research on the 
relationship between peak displacements in the spectrogram and piston errors, the results show a good 
linear relationship. In the close-loop experiment, we make close-loop control for the piston error with 
real-time calculation results and piston error compensating elements, the results show that the system 
remains at the initial interference state with extra perturbation. 

Key words: Piston error detection; Fourier transform; Dispersed fringe sensing; Control 
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